Sensitivity of 7-ray telescopes for detection of magnetic fields in intergalactic medium 
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We explore potential of current and next-generation 7-ray telescopes for the detection of weak 
magnetic fields in the intergalactic medium. We demonstrate that using two complementary tech- 
niques, observation of extended emission around point sources and observation of time delays in 
7-ray flares, one would be able to probe most of the cosmologically and astrophysically interesting 
part of the "magnetic field strength" vs. "correlation length" parameter space. This implies that 
7-ray observations with Fermi and ground-based Cherenkov telescopes will allow to (a) strongly 
constrain theories of the origin of magnetic fields in galaxies and galaxy clusters and (b) discover, 
constrain or rule out the existence of weak primordial magnetic field generated at different stages 
of evolution of the Early Universe. 
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I. INTRODUCTION 

Magnetic fields are known to play an important role in 
the physics of a variety of astrophysical objects, from 
stars to galaxies and galaxy clusters. The existence 
of galactic magnetic fields with strengths in the range 
1 — 10 /iG is established via observations of Faraday ro- 
tation and Zeeman splitting of atomic lines in the radio 
band and of polarization of starlight in the optical band 
P, Q • Magnetic fields of similar strength are found in 
the cores of galaxy clusters Q- Weaker magnetic fields 
with strengths in the range 10~^ — 10~^ G were recently 
discovered at the outskirts of galaxy clusters [1, Q . 

Although the strength and spatial structure of mag- 
netic fields in the Milky Way and some other galaxies 
are reasonably well known today, there is no commonly 
accepted theory about the origin of these magnetic fields 
(see [U, H, 0, 1^ for recent reviews of the subject). There is 
a general agreement that the observed microGauss mag- 
netic fields are the result of amplification of weak " seed" 
fields. The amplification mechanisms under discussion 
are the so-called "a — uj" dynamo (in the case of spiral 
galaxies) and/or compression and turbulent motions of 
plasma during the galaxy/cluster formation processes. 

The nature of the initial weak seed fields for the 
dynamo or turbulent amplification is largely unknown 
P, H, 07 111 • It might be that the seed fields are produced 
during the epoch of galaxy formation by electrical cur- 
rents generated by the plasma experiencing gravitational 
collapse within a proto-galaxy [§, [l^l, or ejected by the 
first supernovae [ll| or active galactic nuclei [l3l- Oth- 
erwise, the seed fields might originate from still earlier 
epochs of the Universe expansion, down to the cosmolog- 
ical phase transitions or inflation times • 

Wide uncertainties in both the mechanism of amplifi- 
cation of the seed fields and in the nature of the seed fields 
themselves have led, over the last half-a-century, to the 



appearance of a long-standing problem of the "origin of 
cosmic magnetic fields" (in galaxies and galaxy clusters) . 

It is clear that the clue for the solution of this problem 
might be given by the measurements of the initial seed 
fields. However, up to recently there was little hope that 
the extremely weak fields outside galaxies and galaxy 
clusters would ever be detected. 

In what follows we show that direct measurements 
of the seed fields and derivation of constraints on their 
nature become possible with the newly available obser- 
vations in the very-high-energy 7-ray band with space 
and ground-based 7-ray telescopes such as Fermi, HESS, 
MAGIC, VERITAS and, in the near future, CTA, AGIS 
and HAWC. The method of measurement of "Extra- 
Galactic" Magnetic Fields (EGMF) with 7-ray tele- 
scopes is based on the possibility of detection of emission 
from electromagnetic cascade initiated by the primary 
7-rays emitted by an extragalactic source and develop- 
ing throughout the InterGalactic Medium (IGM) along 
the line of sight toward the source [H, [H, [isl. [l6l. [l7j. 
Based on the knowledge of sensitivity of existing and fu- 
ture 7-ray telescopes, we find the range of EGMF pa- 
rameters, such as the field strength B and the correla- 
tion length As, in which the EGMF is accessible for the 
measurements with one of the two available measurement 
techniques (imaging [ij, IH, [13 timing [H, [11] of the 
cascade signal). We demonstrate that most of the astro- 
physically and cosmologically interesting range of EGMF 
parameters could be probed with 7-ray observations. 

The plan of the paper is as follows. In Section [IT] we 
summarize the existing bounds on the strength and cor- 
relation length of EGMF which come mostly from radio 
observations. In Section Hill we discuss limits on the cos- 
mological magnetic fields from cosmology. Then, in Sec- 
tion [IV] we compare the existing bounds to the theoreti- 
cal predictions of two classes of models ("astrophysical" 
vs. " cosmological" models) of the "seed" fields and show 
that model predictions normally fall largely below the 
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B, and the correlation length, Ab^. It turns out that 
limits on B imposed by different observations depend on 
A B ■ This means that the limits could be presented as an 
"exclusion plot" in {B,Xb) parameter space, shown in 
Fig. [11 

Magnetic fields in IGM decay due to magnetic diffusion 
over the cosmological time on the distance scales 



^diff = y^~2xf0i3cm (f) 

where Th is the Hubble time and <t ~ fO^^ s~^ is the 
conductivity of the Universe after recombination. This 
means that the correlation length of EGMF is limited 
from below to As > Adiff. At the same time, there are 
no known upper bounds on As and a natural bound is 
only set by the size of the visible part of the Universe, 
A_B < Rh; where Rh is the Hubble radius. The lower 
and upper bounds on Ab arc shown as vertical lines in 

Fig.m 



FIG. 1: Observational limits on EGMF. Cyan shaded region 
shows the upper limit on B imposed by the Zeeman split- 
ting measurement, the lower bound on the correlation length 
imposed by the magnetic diffusion and the upper bound on 
correlation length given by the Hubble radius. Orange shaded 
region shows the limit from Faraday rotation measurements. 
Filled orange region shows the limit derived in the Ref. [2^ . 
while the orange-hatched region is the limit derived in the Ref. 
[2H . Magenta line shows limit which can be imposed by obser- 
vations of deflections of UHECR . Violet vertical-hatched 
regions and the arrows at Ab ^ 0.5 Mpc and As Rh show 
the limits imposed on cosmolo gica lly produced fields by the 
CMB observations [H, HI, liol.liil. Black ellipses show the 
ranges of measured magnetic fields in galaxies and galaxy 
clusters. 



existing bounds. In Sections IVl - Hxl we first summarize 
the methods of measurement of EGMF with 7-ray tele- 
scopes and then estimate the ranges of EGMF param- 
eters which can be probed with different observational 
techniques and different telescopes. Finally, in Section 1X1 
we draw conclusions from our study. 



II. EXISTING LIMITS ON THE EGMF 

Contrary to the magnetic fields in galaxies and galaxy 
clusters, magnetic fields in the IGM have never been de- 
tected. Only upper limits, obtained by different obser- 
vation techniques, exist. In this section we review the 
existing observational limits on the EGMF strength. 

In the simplest settings, the EGMF configuration can 
be characterized by two parameters: the field strength. 



A. Zeeman splitting 

A straightforward upper bound on the EGMF strength 
can be found from the measurements of Zeeman splitting 
of 21 cm absorption line in the spectra of distant quasars, 
which are used to infer the magnetic field in the Milky 
Way galaxy [l^ . The magnetic fields measured via Zee- 
man splitting are usually in the range 1 — 100 ^G and 
are commonly attributed to the field in the Milky Way 
[ist or in other galaxies (see e.g [l^ for detection of 84 
fiG magnetic field in a galaxy at redshift z ~ 0.7). Mea- 
surements of ^ fiG galactic magnetic fields via Zeeman 
splitting technique rule out the possibility of existence of 
still stronger magnetic fields in the IGM. The limit from 
Zeeman splitting measurements, obviously, does not de- 
pend on As and is shown as a horizontal (weakest) upper 
bound on B in Fig. [1] 



B. Faraday rotation 

Measurements of Faraday rotation of polarized radio 
emission from distant quasars provide a possibility of de- 
tection of EGMF of the strength somewhat lower than 
the one accessible for the " direct" measurements via Zee- 
man splitting. The rotation measure RM = Ax/AA^ 
(Ax is the change of the polarization angle between the 
wavelengths A and A -I- A A) is proportional to the product 
of magnetic field component along the line of sight, B^ 



A third important parameter is the volume fiUing factor V of 
magnetic fields of a given strength B and correlation length 
Omitting this parameter we restrict ourself to the task of the 
search of the "dominant" EGMF, with volume filling factor V ~ 1 
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and the electron density of the IGM rig 



Deflections of UHECR 



Ax 
AA2 



Q.2°T{ze) 



1,0 



10-11 G 



«e,0 



10-7 cm-3 



(2) 



where z^; is redshift of the source of polarized emission 
(we use index E for "emitter") and function T{ze) — 
Ho J^"^ {1 + z)^ {dt / dz)dz , if one assumes a simple redshift 
evolution of B^^{z) ^ i?||.o(l + nf.{z) = ne,o(l + zY 



by 



In the ACDM cosmology, the derivative dt/dz is given 



dt 
dz 



Hail + z) ^nM{l + z)3 + QA 



(3) 



with Ho ~ 71 km/s/Mpc, Qm ^ 0.27 and Qa ^ 0.73 be- 
ing the present-day value of Hubble parameter, the cold 
dark matter and cosmological constant energy densities, 
respectively [20| . 

Under simplest assumptions, free electrons are 
distributed homogeneously through the intergalactic 
medium. Using this model for electron distribution, a 
limit of B < 10-11 G was obtained in the Ref. [2l| for 
the present-day electron density close to the critical den- 
sity of the Universe. Re-scaling this limit for the electron 
density of the order of the baryon density, one would find 
a limit which is a right-bottom corner of the hatched or- 
ange region in Fig. [TJ 

The limit on EGMF imposed by the Faraday rotation 
measurements depends on the EGMF correlation length 
Xb- If As <C Rh, the polarization angle experiences 
random changes due to the passages of multiple "cells" 
of the size R ^ Xb with coherent magnetic field. This 
means that the x changes proportionally to the square 
root of the distance. 



Ax 
AA2 



Vdx'i 



(4) 



where D is the distance to the radio source. Non- 
detection of the EGMF induced Faraday rotation im- 
plies, therefore, an upper limits on B which scales as 



B 



Xg^^. This upper limit is shown in Fig. [T] as an 
orange-hatched region. 

The simplest assumptions about homogeneous distri- 
bution of free electrons in the IGM might be an over- 
simplification. More complicated models of electron dis- 
tribution were considered in the Refs. [1^ and [1^. In 
these references IGM models based on the Lya data were 
considered. Under certain assumptions about the depen- 
dence of Tie and B on the density of the Lya clouds, the 
authors of Ref. [1^ derive a limit on B which is weaker 
than the one cited above (shown as a filled orange region 
in Fig. [1]). The limit on B\\ depends weakly on Xb if 
Xb 3> Aj, where Aj is the Jeans length scale, which char- 
acterizes typical distance between the Lya clouds [1^ . At 
small correlation lengths, Xb ^ Aj, the limit scales as 

— 1/2 

B|| ~ A^ since the Faraday rotation angle experience 
random changes during the passage through each cloud. 



Magnetic fields in IGM could be probed by the mea- 
surements of their effect on trajectories of charged par- 
ticles (high-energy electrons and cosmic ray protons or 
nuclei), if their sources are known. Defiections of high- 
energy electrons and positrons are discussed below in 
Section |Vl In principle, measurements of deflections of 
Ultra-High Energy Cosmic Rays (UHECR) with ener- 
gies i^uHECR > 10^^ eV also can be used to constrain 
the intergalactic magnetic fields. In the simplest case, 
when the energy losses of UHECR could be neglected, 
the defiection angle of UHECR in a regular magnetic 
field with coherence length Ar larger than the distance 
to the source D is given by [2J| : 
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where where B±^ is the magnetic field component orthog- 
onal to the line of sight, £'uhecr is particle energy and 
Z is the atomic charge. 

In an opposite case Xb ^ D one has [23 |: 
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The main energy loss channels of UHECR protons is 
pion production on CMB [2^ , while for heavy nuclei it is 
photo-disintegration on cosmic infrared background [2^ . 
In both cases the energy / charge attenuation distances 
arc < 100 Mpc. This limits the distances toward the 
sources of the highest energy cosmic rays to be not larger 
than D - 100 Mpc. 

It is clear that measurement of EGMF with UHECR 
would be possible only under the condition that ex- 
tragalactic point sources of UHECR would be detected 
(which is not the case at present). Even if the cxtragalac- 
tic UHECR point sources would be known, attempts of 
the measurement of EGMF would face the following prin- 
cipal difficulty. For a known UHECR source, deflections 
of UHECR arrival directions from the real source posi- 
tion are determined not only by EGMF, but also by the 
defiections in the Galaxy, by magnetic fields in the in- 
tervening large scale structure elements, like galaxies or 
galaxy clusters along the line of sight and by deflections 
in the source host object (galaxy, galaxy cluster). As a 
result, the defiection angle is a sum of at least three 
terms. 



^Gal + SeGMF + ^'source 



(7) 



Measurement of deflections by EGMF, 6'egmf, via mea- 
surement of 9 implies that the deviation by the Galactic 
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magnetic field. 6'Gai and possible deviations by the source 
host or intervening galaxies / galaxy cluster(s), ^sourccj 
arc known. 

The magnetic field of the Milky Way galaxy is con- 
ventionally modeled as a sum of the regular and turbu- 
lent components of the field in the disk and halo of the 
Galaxy. This means that the term ^cai in the above de- 
composition of 9 is itself a superposition of at least four 
terms: 



JGal 



qrcgular 
'^Disk 



qturbulciit 
^Disk 



qrcgular 
-^Halo 



qturbulciit 
^Halo 



(8) 



Deflection 
of Galactic 



by regular and 
Disk and Halo, 



turbulent components 

/iDisk /iHalo ct-nA 

'^regular' '^regular 

^?urbuicnt' ^'uirbuicnt ' Can be estimated by substitution of 
the typical disk/halo size at the place of D and of the typ- 
ical disk/halo field strength and correlation lengths at the 
place of B,Xb in Eqs. (0), dH). Deflections of UHECR 
by the regular field in the disk were studied in 

many theoretical models staring from Ref. [l^]. Typi- 
cal values of parameters entering the analog of Eq. ([5|) 
imply substitution D — > 2 kpc length scale (for sources 
located far from the Galactic Plane) and B ^ 2 fiG, 
which give for Eq. dS]) 6*0?;;'^' ~ 4°. Turbulent fields are 
typically assumed to have coherence scale ~ 50 pc and 
B ~ 4 ^G, which, for same field height scale ~ 2 kpc 
gives 6lj^7^5j"i™* ~ 0.5°. Contributions of the Halo fields 
are less certain, but result in deflections at least of the 
same order, see recent discussion of all components in 

Refs. Illia,!!^. 

Although deflection angles of UHECR by the Galac- 
tic magnetic field can be readily estimated by order of 
magnitude, uncertainties in the measurements of Galac- 
tic magnetic field and discrepancies between the existing 
measurements and existing theoretical models of Galac- 
tic fields [2^ H^l do not allow to predict the deflection 
angle and direction of deflection for particular lines of 
sight (toward UHECR sources). This means that, most 
probably, the details of the structure of the Galactic mag- 
netic field along the line of sight toward UHECR sources 
would have to be deduced from the UHECR data itself, 
rather than just taken into account in the UHECR data 
analysis [3l| . This, obviously, will introduce a large un- 
certainty into the derivation of the properties of EGMF 
from the UHECR data. 

Finally, the last term in the Eq. ([7]), 6'source, is equally 
uncertain. Recent attempts of modeling of deflections of 
UHECR by the intervening elements of large scale struc- 
ture, such as galaxy clusters and/or filaments by two 
groups (sec Refs. [s^ and [s^) give contradictory re- 
sults, which reflects uncertainties of the structure of mag- 
netic fields inside and around clusters and filaments. In 
addition, if an UHECR source is nearby, the source host 
galaxy or galaxy cluster could span several degrees on the 
sky. Significant deflections of UHECR by magnetic fields 
in the host galaxy or galaxy cluster could produce intrin- 
sic 1-10 degree-scale extension of the source [U, which 
would make the extraction of information about EGMF 



from the study of deflections of UHECR arrival directions 
from the source position still more problematic. 

Neglecting the above-mentioned problems, one could 
derive a theoretical "sensitivity" limit of future UHECR 
experiments for the detection of EGMF. Taking into 
account the typical angular resolution of current and 
next generation UHECR experiments, like JEM-EUSO 
is 6'psF ~ 2° [H, one could find, from Eqs. © and © 
that EGMF with the strength down to S - 10"^" G 
(in the case of large correlation length As) could infiu- 
ence the observational appearance of the signal from the 
UHECR source. Theoretically possible sensitivity limit 
of UHECR experiments for the measurements of EGMF 
is shown as magenta line in Fig. [T] 



III. LIMITS ON THE COSMOLOGICAL 
MAGNETIC FIELDS FROM COSMOLOGY 

Zeeman splitting or Faraday rotation methods allow to 
detect, or put upper bounds, on the weak "seed " mag- 
netic fields in the intergalactic medium in the present 
day Universe. As it is mentioned in Introduction, such 
fields could have been generated either during the epoch 
of galaxy formation, or at the earlier stages of evolution 
of the Universe. If the seed fields were generated be- 
fore the epoch of recombination, additional constraints 
on the field strength and correlation length can be de- 
rived from the analysis of the Cosmic Microwave Back- 
ground (CMB) data. 

Most of the cosmological " magnetogenesis" models re- 
sult in predictions of tangled magnetic field configura- 
tions with broad band power-law like spectra in Fourier 
space, cut off at a characteristic (model dependent) 
length scale As, so that the Fourier components of the 
field have the form [30| 



Air 



-e{kB-k) 



(9) 



where ks ~ 2tt/Xb and n is the power-law index. The 
energy density contained in the EGMF in the present day 
Universe, — Bq/{8tt), is expressed as an integral over 
the Fourier space p% = l/(87rfc|) / d^k\Bk\'^. Note that, 
in order to have finite energy in magnetic field, one has 
to assume n > — 3 in the Eq. 

Constraints on the parameters of the magnetic 
field spectrum come from the measurements of the 
anisotropics of the CMB spectrum, Faraday rotation of 
CMB, chemical potential of CMB and from Compton pa- 
rameter measurements (analog of Sunyaev-Zel'dovich ef- 
fect at recombination epoch) . If the magnetic fields were 
generated prior to the epoch of Big Bang Nucleosynthesis 
(BBN), additional (weaker) constraints could be derived 
from the BBN calculations, at the level of Bq < 10^^ G 
(see, for example, Q)- It is important to note that con- 
straints on cosmological magnetic field strength B nor- 
mally depend not only on the correlation length A b , but 
also on unknown power law index n. 
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The most straightforward limit on the strength mag- 
netic fields homogeneous over the Hubble distance scale 
could be derived from the non-observation of the large 
angular scale anisotropics of the CMB. This gives an up- 
PCT limit B < 4x 10"^ G for the fields with Xb ^ Rh 
[37| . The powerlaw index dependent limit on the mag- 
netic field strength were derived from the analysis of the 
CMB angular power spectrum in the Refs. 38, 39, .40j . 
The envelope of the upper bounds on {B,Xb) for the 
range of the powerlaw indices — 3 < n < 2 is the lower 
boundary of the thin vertical violet hatched region of 
{B, Xb) parameter space marked "CMB" in Fig. [1] 

Non-thermal dissipation of magnetic field energy into 
the energy of electrons/positrons during the recombina- 
tion epoch could lead to distortion of blackbody CMB 
spectrum [4l| . This distortion produces non-zero chemi- 
cal potential /i, which was calculated in the Ref. [36] in 
the form of a double integral. In order to calculate the 
bounds on B, Xb, we consider two limiting cases in which 
the integral could be taken analytically, namely, the 
cases when the magnetic field correlation length is much 
smaller or much larger than the characteristic damping 
length scale (the wavelength at which the magnetic field 
is damped by a factor e): 



An - 



to 



2n 



3/2 V 15n"aT/. 



400 pc. 



(10) 



where time constant is to ~ 2.4 x 10^^ s, is the charac- 
teristic redshift of freeze-out from double-Compton scat- 



tering Zfj, — 2.5 X 10° 



is electron density and axh is 



Thomson cross section. 

Taking into account the constraint < 9 x 10~^ at 95 



% confidence level from COBE FIRAS data [43, one can 
derive limits on magnetic fields in the analytical form in 
the two limiting cases. In the case ^ Xo one has [1^: 



B < 3.2 X lO^^'G 



1 



A, 



K V400 pc 



-(n+3)/2 



(11) 



where K = 1.4r(7i/2 + 5/2)r(3n/5 + 
9/5)2-("+5)/2(6/5)(7i + 3) is a constant of the or- 
der of unity, K = 0.8/2.1 for n = -2/ + 1. In the 
opposite case As ^ A_d the constraint is 



S < 3.2 X lO^^G 



1 



A, 



400 pc 



(12) 



where K2 is another constant of the order of unity. The 
strongest limit on the field strength is obtained for the 
fields with correlation length Xb ^ Xd- For the fields 
with such correlation length no convenient analytical ap- 
proximation could be found and instead a numerical inte- 
gration of expression of the Ref. [11] should be performed. 
In Fig. [1] we show the bound on B,Xb implied by the 
analysis of distortions of CMB spectrum as extrapolation 
of the analytical approximations given by Eq. pT|) . (fT^ 
for the whole ranges A^ < A^i, Xb > Xd- Note that the 



dependence of the limits on B,Xb on the power-law in- 
dex n practically disappears in the case Xb ^ Xd (only 
weak dependence remains in the constant )■ At the 
same time, uncertainty in the value of n "washes out" 
the upper bounds on B at the length scales Xb ^ Xd- 

We have cross-checked the above results by adopting 
an assumption that the magnetic field spectrum has the 
(5- function rather than a power-law shape of Eq.([5|). In 
terms of Eq. Q this would correspond to the limit n — > 
00. As expected, Eq. (jll[) does not give any constraints 
in this case for Xb < Xd, while limit Eq. (|12p remains as 
it is. 

Apart from producing a non-zero chemical poten- 
tial, transfer of the magnetic field energy to electrons/ 
positrons could result in non-zero Compton parameter 
y. Taking into account restrictions y < 1.5 x 10~^ from 
COBE FIRAS^ne finds a limit B <3x 10'^ G at A - 
0.3-0.6 Mpc [si]. Note, that limits Eq. ^ and Eq. ^ 
constrain magnetic fields created at z > = 2.5 x 10^ 
, while the limit following from restrictions on y applies 
for fields created before recombination z > 2 x lO"'. 

Magnetic fields created before recombination could 
produce another observable phenomenon: Faraday ro- 
tation of linear polarization of CMBR. Constraints on 
B, Xb stemming from non-observation of this effect were 
first discussed in the Ref. and subsequently updated 
using the 5-years data of WMAP in the Ref. [ij]. The 
limits coming from non-observation of Faraday rotation 
in CMB signal are, at present, weaker than the limits im- 
posed by the rotation measures of distant blazars or lim- 
its from the CMB angular power spectrum [i^] . The limit 
from non-observation of Faraday rotation in the CMB 
signal is shown by thick violet vertically hatched region. 



IV. THEORETICAL PREDICTIONS 

Search for extremely weak EGMF is important in the 
context of the problem of the origin of magnetic fields in 
galaxies and galaxy clusters. Most commonly accepted 
hypothesis about the origin of magnetic fields in spiral 
galaxies is that they were produced via the so-called 
"a — w dynamo" mechanism from tiny "seed" fields of 
uncertain origin [l| . Estimates of the efficiency oi a — lu 
dynamos imply that the " seed" field should have strength 
of the order of ~ 10^^^ G or higher, although these es- 
timates suffer from large uncertainties^. Alternatively, 
magnetic fields in galaxies and galaxy clusters could be 
produced via compression of the seed magnetic fields 
present during structure formation epoch, and their am- 
plification by the turbulence (see e.g. Ref. [3). Decisive 
observational test of alternative theories of the origin of 



^ Recent discoveries of strong magnetic fields in galaxies at signifi- 
cant rodshifts imply still stronger "seed" magnetic fields, see e.g. 



El 



6 




-12 11- 10-9 -8 -7 -6 -5 -4 -3 -2 1 1 2 3 4 12 11 10 9 B 7 6 5 4 3 2 1 1 2 3 4 -12-11-10-9 -8 -7 -6 -5 -4 -3 -2 -1 1 2 3 4 
log(A, [Mpo]) log{K [Mpo]) log(A, [Mpo]) 



FIG. 2: Model predictions and estimates for the EGMF strength. Cyan shaded region and black elhpses show the experimental 
limits and measurements from Fig. [5] Upper bound at B ~ 10"^" G shown by solid line comes from flux conservation during 
galaxy formation argument J^^- Upper bound at B = lO"^'^ G shows a limit imposed by constrained simulations of magnetic 
fields in galaxy clusters [l^. |33]. Left panel: left and right hatched regions show theoretically allowed range of values of B, Xb 
for non-helical and helical fields generated at the epoch of electroweak phase transition during radiation-dominated era. Middle 
panel: left and right hatched region show ranges of possible B, Xb for nonhelical and helical magnetic fields produced during the 
QCD phase transition. Right panel: hatched region is the range of possible B,Xb for EGMF generated during recombination 
epoch. Dark grey shaded region shows the range of {B,Xb) parameter space accessible for the 7-ray measurements via 7-ray 
observations. Light-grey shaded regions show the parts of the parameter space in which the existence of EGMF could be 
confirmed or ruled out, but no measurements of EGMF strength is possible. 



cosmic magnetic fields can be given only by observations 
of the initial weak seed fields implied in all the theoretical 
models. 

However, the nature of the weak seed fields remains 
largely unconstrained. Theoretically, two main hypothe- 
ses exist: that of the " astrophysical origin" (see e.g. [l|) 
and of the " cosmological origin" (see e.g. 0]) of the seed 
fields. The " astrophysical origin" models usually assume 
that the seed fields are produced via the so-called " Bier- 
mann battery" effect (difference in mobility of electrons 
and protons in astrophysical plasma) at the early stages 
of galaxy formation. Different types of collective motions 
of plasma could be involved, such as the outflows from 
the first supernovae activity of active galactic nu- 
clei 12 1, gravitational collapse 0, [l^ and/or turbulence 
[H; 47 1 in proto-galaxies. The "cosmological" models 
usually exploit a similar effect at earlier stages of evo- 
lution of the universe, during or before recombination 
m, m, electroweak 0, and QCD [M H p| 

phasetransitions or at still earlier epochs [sgI. ISTI-ISsI. 59 j 



(see Q for a review). It is clear that detection of the 
seed magnetic fields would not only help to resolve the 
problem of the origin of magnetic fields in galaxies and 
galaxy clusters, but also provide a new observational data 
constraining physical conditions in the Early Universe. 

Significant difference between the "cosmological ori- 
gin" and " astrophysical origin" models of the seed fields 
is that in the former case weak seed magnetic field should 
be present everywhere the Universe today. In particu- 
lar, the seed fields could be found in the voids of large 
scale structure, outside galaxies and galaxy clusters. To 



the contrary, if the seed fields are produced via astro- 
physical mechanisms inside proto-galaxies, no magnetic 
field generation outside galaxies and clusters is expected. 
Magnetic field in the voids of the large scale structure is 
expected to be close to zero. 

The "cosmological origin" models could be divided 
onto two broad classes: models based on mechanisms 
operating during phase transitions in the Early Universe 
and models based on mechanisms operating during infia- 
tion epoch. In the latter case, the fields might be pro- 
duced on a super-horizon scale and no firm theoretical 
limits on the characteristic field correlation length can 
be deduced. In the former case, the strength and corre- 
lation length of cosmologically generated magnetic fields 
are well limited from above. 

The correlation length of magnetic field can not exceed 
the size of the horizon at the moment of phase transition 



A, 



100 GeV 



cm (13) 



where a, H and T are scale factor, expansion rate and 
temperature of the Universe and Gn is the Newton con- 
stant. Indices and * refer to the parameter values to- 
day and at the moment of magnetic field production, 
respectively. The magnetic field strength at the corre- 
lation length As is limited by the requirement that the 
magnetic field energy density pB = B"^ /(%'k) should not 
overdose the Universe. Since the magnetic field energy 
density evolves in the same way as the radiation energy 



density, pB ^ Pph 



present day magnetic field at 
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the scale Xb is limited to be below the field for which the 
energy density is equal to the CMB energy density, 

B{Xb,) < V47rg*PcMB ^ 3 X 10-«(.g,/2)i/2 G (14) 

where is the number of relativistic degrees of freedom 
at the moment of magnetogenesis. If the magnetic field 
would grow so that it would start to contribute signifi- 
cantly to the overall energy density of the Universe, its 
further amplification during the magnetogenesis epoch 
would stop because the energy transferred to the mag- 
netic field would be dissipated via strong gravitational 
wave production [iO, [l^l • 

The upper limit on the magnetic field strength at the 
horizon scale (jl4p implies an upper limit on super-horizon 
scales which follow from straightforward causality argu- 
ments [6l| . For non- helical magnetic fields an upper limit 
on the strength of magnetic field at arbitrary length scale 
A could be obtained taking into account the fact that the 
power-law index of magnetic fields produced in a causal 
way at phase transition(s) arc limited to be n > 2 [61 1 



B{Xb) < BiXB.) ( ^ 

^B* 



h3)/2 . , s -5/2 



< 3x10" 



As 
Xb* 



(15) 

Energy contained in magnetic fields with small enough 
correlation length is dissipated in the course of evolution 
of the universe. This leads to the increase of the "integral 
scale" (i.e. the distance scale which gives dominant con- 
tribution to the magnetic magnetic field energy density) 
with time [6^ up to the scale 



A 



B,I 



V 

H 



(16) 



where v is the characteristic velocity scale, which is of 
the order of either Alfven or viscous velocity at dif- 
ferent epochs of the Universe expansion. Numerically, 
Xb,i ~ 1 [-B/5 X 10~^^ G] kpc for magnetic fields pro- 
duced much before recombination epoch and Xbj ^ 
1 [-8/8 X 10^^^ G] kpc for the magnetic fields produced 
at recombination. 

The above constraints limit the possible present day 
values of B,Xb for the magnetic fields produced in the 
Early Universe. The possible ranges of the {B,Xb) pa- 
rameter space are shown in Fig. [2] for the cases when 
magnetogenesis proceeds during electroweak or QCD 
phase transitions or at the moment of recombination. 
These regions are bound on the left by the constraint on 
Xbj given by Eqs. ([TB)) . The right side boundaries 

of the allowed regions are determined by the constraints 

(HI, (HH). 

If the magnetic fields possess non-zero helicity, the en- 
ergy contained in the short-wavelength modes can be 
transferred to the longer-wavelength modes via the de- 
velopment of "inverse cascade", rather than dissipated 
in the fluid motions. This results in a slower decay of the 
magnetic field energy density and faster growth of the 
integral scale, than in the case of non-helical fields [6^ . 



For the helical magnetic fields, the magnetic field 
power can be transferred from the initial scale Xb* to 
a larger scales Xb via inverse cascade. Conservation of 
helicity implies that H ^ B^L, where L is the distance 
scale, is conserved during the inverse cascade. Analyt- 
ical/numerical calculations of the inverse cascades show 
that the cascade Icngthscale evolves as aL ~ t^/'^ [g^l 
during the radiation dominated epoch, which means that 



A, 



B{Xi 



_T*_ 
1 eV 

1 eV 



1/3 



A, 



-1/6 



B{Xi 



(17) 



at the end of the radiation-dominated epoch at T ~ 1 eV. 
The causality argument applied for the helical magnetic 
fields gives the constraint on the power-law index n > 3 
at the scales A > As [HH. The limit (fTT]) is weaker than 
the bound (|15p on non-helical magnetic fields. This limit 
determined the right-side boundary of the allowed regions 
of {B, Xb) parameter space for the electroweak and QCD 
phase transition magnetogenesis in Fig. [2l 

It is interesting to note that predictions for the 
strength and correlation length of the " primordial" mag- 
netic fields produced at electroweak and QCD phase tran- 
sitions fall in a region of {B, Xb) parameter space which 
is not accessible for the existing measurement techniques, 
such as Faraday rotation or Zeeman splitting methods. 
However, it turns out that this region of {B, Xb) param- 
eter space is accessible for the measurement techniques 
which exploit the potential of the newly opened field of 
very-high-energy (VHE) 7-ray astronomy. In the fol- 
lowing sections we demonstrate that using the current 
and next generation ground and space based 7-ray tele- 
scopes one could probe the part of (B, As) parameter 
space shown by the light and dark grey-shaded regions 
in Fig. [5] and in this way test the models of the origin of 
seed fields and the models of the origin of magnetic fields 
in galaxies and galaxy clusters. 



MEASUREMENTS OF EGMF WITH 7-RAY 
TELESCOPES 



A. Basic formulae 

Multi-TeV 7-rays emitted by distant point sources are 
not able to propagate over large distances because of 
the absorption in interactions with optical/infrared ex- 
tragalactic background light (EBL). 

Redshift-dependent inverse mean free path of such 
gamma-rays of the energy E'^^ propagating at redshift 
z through the EBL of the density UEBbi^, z) can be de- 
fined as 



de' 



< (T-y-,nEBL > 



dn 



EBL 



de' 



dfi{l - ti)ay.y{s), (18) 
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where (J^j{s) is pair production cross section which de- 
pends on s = 2E'^^e'{l — /x) with /i being cosine of 
the angle between the directions of propagation of 7- 
ray and background photon. The Umit of integration 
^'min = ''^e/^'-yo (O corrcsponds to the threshold 

of the pair production. 

Estimates of D^[E'^^,z) suffer from significant un- 
certainty in the modeling of nEBL{(-,z). There ex- 
ist several theoretical models for the formation of in- 
frared/optical background based on the models of evo- 
lution of starlight and dust emission from different types 
of galaxies at different redshifts. The resulting estimates 
oi UEBLi^, z) widely differ both in the predicted shape of 
local UEBLi^, z = 0) and in its evolution with increasing 
redshift, sec Refs. [gI, [b^, IH, . However, all the mod- 
els agree in a general trend of decrease of with the in- 
creasing energy and redshift. Following results of ref. [67j . 
one can assume that UEBhi^, z) « (1 + z)^'^nEBL{^, z = 
0). In this case the integrals in the Eq. (|18p depend on z 
only through E'^^^ . Simplifying the expression in this 
way, we have found that in the energy band of interest 
for the following discussion, a broad range of model pre- 
dictions/uncertainties for Dj could be described by the 
following simple approximation 



i?^(i?;„,z) = 4o 



[i + zy 



20 TeV 



Mpc 



(19) 



where a numerical factor k = 
for the model uncertainties, 
calculations of the Refs. [13 
that in the energy range E'^^ 



= K{Ejg,z) ~ 1 accounts 
Comparing the results of 
[6^ . ISSl . [66j one can find 
~ 0.1 - 10 TeV and the 



redshift range z < 1 the uncertainty in k is as large as 
0.3 < K < 3. 

As the 7-ray propagates from the source toward the 
Earth, the optical depth with respect to the pair produc- 
tion grows as 



dr _ 1 

'dt ^ D^{E' z) 



(20) 



r reaches 1 at the time <^-y, which can be implicitly found 
from equation 



dt 



D,{E'^„.z) 



dz 



dt 



D^{E'^^^,z) dz 



(21) 



where is the time of emission of photon from the source 
and ze and z^^ are, the redshifts corresponding to the 
times Ie and t^^. 

In principle, the mean free path of very-high-cnergy 7- 
rays from distant sources could be large enough so that 
z^^ could be significantly different from ze- In this case 
the exact expression pip should be used to determine 
t^^. The "comoving" mean free path d^ of 7-rays of 
initial energy E' (ze) could then be estimated as 



ait') 



(22) 




FIG. 3: Geometry of propagation of the direct and cascade 
7-rays from the source (on the left) to the observer (on the 
right). 



where a(t) = ao/(l + z). 

In the following calculations (which are mostly order- 
of magnitude estimates) we will adopt a simplifying as- 



sumption that z. 



77 



Ze- In this case t. 



be explicitly found from Eq. 
reduces to 

d, [E'^^izE)] 



tE could 
2T|) and expression ((22l) 



77 



D-^jE'^^.ZE) 
a{tE) 



(23) 



One should remember, however, that assumption z-^^ ~ 
Ze does not hold for relatively low energy (sub- TeV) 7- 
rays which could produce cascade signal in Fermi (GeV) 
energy band. 

Using Eq. ^ one can find the optical depth with 
respect to the pair production for the 7-rays of the energy 



that reach observer on the Earth: 



t{E^„,ze) 



, dt 

d. 



dz D^{{1 + z)E^„,z) 



(24) 



Adopting the approximation p9p for one finds the 
estimate for r 



t{E^oiZe) 



2(v/r>M(i + z)3 + r>A-i) 

3HoCtMD^{E^„0) 



(25) 



Uncertainties in the model predictions for D^{E'^^^^ z) re- 
sult in the discrepancies between the model predictions 
for r. Comparing the results of different calculations re- 
ported in the literature, one can find, for example, r = 1 
at z = 0.03 for E^ = b-l% TeV in the models of rcf. [H, 
= 9 TeV in the model of je!] , -B^ = 7 - 8 GeV in the 
models of ref. 16511 and E^ = 2.7 TeV for fast evolution 
model of ref. [6y|, the value from Eq. (PS)) is i?^ = 6 
TeV. One can verify that the estimate of r given by Eq. 
^5l) is within factor 2 from results of [64| for 2 < 1 and 



0.2 TeV < < 30 TcV. Taking into account the fact 
that error introduced by approximation (j25p is smaller 
than uncertainty of the models, we will use approxima- 
tions of Eqs. and Eq. (pS)) instead of exact Eqs. (fT^ 
and Eq. ([24]) in the following sections. 

The pair production on EBL reduces the flux of 7-rays 
from the source by 



F{E,„)^F,{E'{zE))e-^^'^- 



(26) 



where F{E^g) is the detected spectrum, Fq{E'^^) is the 
initial spectrum of the source and t^Ej^, ze) is the opti- 
cal depth (jM)) . (pS)) . 
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The e~^e~ pairs produced in interactions of multi- 
TeV 7-rays with EBL photons produce secondary 7-rays 
via inverse Compton (IC) scattering of the Cosmic Mi- 
crowave Background (CMB) photons. Typical energies 
of the IC photons reaching the Earth are 



^ (1 + Z^jj 



E 

) ^CMB~~2 — ^■'^'^ 



EL 



TeV 



20 TeV 

(27) 

where e'c^./^ = 6 x 10 "'(I + z^^) eV is the typical en- 
ergy of CMB photons. In the above equation we have 
assumed that the energy of primary 7-ray is E'^^^ ~ 2E'^ 
with E'^^ being the energy of the primary 7-rays at the 
redshift of the pair production. Upscattering of the in- 
frared/optical background photons gives sub-dominant 
contribution to the IC scattering spectrum because the 
energy density of CMB is much higher than the density 
of the infrared/optical background. 

Deflections of e+e~ pairs produced by the 7-rays, 
which were initially emitted slightly away from the ob- 
server, could lead to "redirection" of the secondary cas- 
cade photons toward the observer. This effect leads to 
the appearance of two potentially observable effects: ex- 
tended emission around an initially point source of 7-rays 
[1^ . [l6l . [17} and delayed " echo" of 7-ray flares of extra- 
galactic sources [l^ [l5| . 

The cascade electrons loose their energy via IC scat- 
tering of the CMB photons within the distance 



10^^(1 + z^^) 



El 



The deflection angle of the e^e 



10 TeV 

(28) 

pairs, accumulated over 



the cooling distance, depends on the correlation length 
of magnetic field, \b- 

If the correlation length As is much larger than _De, 
motion of electrons and positrons at the length scale Df, 
could be approximated by the motion in homogeneous 
magnetic field. In this case typical deflection angle 5 
is estimated as a ratio of to the Larmor radius in 
magnetic field B' , 



Rl 



eB' 



3x10 



28 



" B' ■ 


-1 


\ 1 


10-18 G 




10 TeV 



cm. (29) 



Note, that, in principle, EGMF depends on the redshift, 
B' = B'{z). In the simplest case, when the magnetic 
field strength changes only in result of expansion of the 
Universe, B'{z) ^ i?o(l -I- z)^, where Bq is the present 
epoch EGMF strength. This gives 



-De 

Rl 



3 X 10~^(H-z-y^)" 



B' 



3 X lO-^\\ + z^^) 



Bn 



10-18 G 
El 



EL 



10-18 G 



10 TeV 

-2 



10 TcV 



(30) 



If the correlation length of magnetic field, , is much 
less than electron cooling distance D^, electron deflec- 



tions are described as diffusion in angle, so that the de- 
flection angle is estimated as 



S = 



Rl 
B' 



5 X 10-^(1 + 2.^.^) 



EL 



10-18 G 



1 kpc 



1/2 



5 X IQ-^l + z^^) 



-1/2 



EL 



10 TeV 



-1 -3/2 



-3/2 



(31) 



10 TeV 



10-18 G 



-,1/2 



1 kpc 



where we have assumed that A^ scales with 2; as A^ = 
Aso(l + z)~^ with Abo being the present epoch EGMF 
correlation length. 

Knowing the deflection angle of electrons, one can 
readily find the angular extension of the secondary IC 
emission from the e+e~ pairs using simple geometrical 
calculation in the comoving reference system, shown in 
Fig. [3] In this figure is the comoving distance to the 
source, defined as 



ao-ffo Jo V^m(1 + zf + flA 



zdz 



(32) 



*^ dt 

and dj is the 7-ray mean free path given by Eq. (j22|) 
which, in the case z^^ ~ ze can be calculated using Eq. 
(f23ll . The angle 9ext is expressed through dE,d^ and 6 
as 



sin(eext) 



sin (5 



(33) 



to 



In the the above expression reduces 

D-,iEL,„,z)^ S 



De{z) 



-5 = 



T0{Ej^,z) 



(34) 



where Dg = a{tE)dE is the angular diameter distance, 
Tg = Dq/D^ and we have assumed tq > 1,(5 ^ I. For 
small z, Tg is close to T{Ejg,z). Numerically, the above 
estimate of the size of extended emission around an ex- 
tragalactic point source is 



0.5°(1 

E^ 



0.1 TeV 



10 J 



Bn 



10 



-14 



G 



o.o7°(i + z)-i/2 [i^; 





-3/4 


r Bn 1 




Abo 


0.1 TeV 




Lio^i* gJ 




1 kpc. 



1/2 



A'r > De 



(35) 

A clear observational signature of the presence of EGMF 
induced cascade emission around an initially point source 
is the decrease of the extension of the source with the 
increase of photon energy. 
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Difference in the path between the direct and cascade 
7-ray leads to the appearance of delayed emission from 
a flaring extragalactic 7-ray source. The geometrical 
scheme of Fig. [3] enables to estimate the time delay be- 
tween the direct photons and secondary inverse Compton 
7-rays produced by deflected electrons as 



^dclay — ^77 {^E ^cascade) 



(36) 



where tj^ is found from Eq. (|2T|) . tE is the light travel 
time of direct 7-rays from the source, 



tE 



dt , 
—az 

dz 



(37) 



and ^cascade is the travel time of the secondary IC pho- 
tons, which is implicitly expressed through the equation 



dt' 



{d 



di 



2dEd^^ cos((5 + Ocxt)) 



1/2 



where to is the present time 

In the case z. 
sion is 



(38) 



T, 



delay 



(1 + ^) 



ZE = z, 6 <^ 1 the resulting expres- 

D,{E'^„,z)S' j,(j;;,z) \ 

De{z) J 



Numerically, at z ^ 1, the time delays for the cases of 
magnetic fields with large and small correlation lengths 
are given by 



T, 



delay 



7 X 10'^k(1 



V^)(l + z)- 



E-y 


-5/2 


Bo 


0.1 TeV 




Uo-18 gJ 



s, 



lO^K{l~Tg'){l + Z)- 





-2 


\ ^" 1 


2 


Abo 


0.1 TeV 




Lio-^« gJ 




1 kpc 



X'n > De 



From the above equation one can see that lowering the 
energy threshold of a 7-ray telescope enables to measure 
weaker magnetic fields. This means that the weakest 
magnetic fields at the level of the lower bound for the 
galactic dynamos, B ~ 10~^^ eV could be best explored 
with Fermi which detects photons with energies above 
0.1 GeV. 



VI. LOWER LIMIT ON EGMF STRENGTH VS. 
"UNIVERSAL GEV TAILS OF TEV 7-RAY 
FLARES" EFFECT 

From the discussion of Section IIVI it is clear that a 
possibility that the EGMF strength is B ^ -Biimit, where 
^limit is given by Eq. (|4ip . is not ruled out observation- 
ally and/or theoretically. Moreover, most of the "astro- 
physical origin" models of the seed fields for the dynamo 
/ compression / turbulence amplification invoke mech- 
anisms which produce seed fields locally in the parent 
(proto)galaxy, rather than throughout the IGM, so that 
the field in the IGM is expected to be nearly zero. 

The possibility of the absence of magnetic fields in ex- 
cess of -Biimit (|4T|) in the IGM could be readily tested via 
observations with Fermi. Indeed, if the EGMF in the 
voids of large scale structure does not significantly de- 
flect e"'"e~ pairs produced via absorption of the highest 
energy 7-rays from TeV blazars, one still expects that 
the IC emission from the secondary e"*"e~ pairs will pro- 
duce an observable contribution to the detected 7-ray 
signal. Moreover, the cascade emission is expected to 
produce a delayed "tail" of the 7-ray flares, because of 
the deflections of e+e^ pairs by an angle ^umit ~ rne/E,,. 
Substituting (Sumit at the place of S into Eq. ([5^ one 
finds 



.(40) 

The magnetic field induced time delays could be iden- 
tified in the observational data via their characteristic 
dependence on the photon energy. 

Magnetic field induced extended emission and delayed 
" echo" of 7-ray flares are detectable only if the deflection 
angle S is larger than the intrinsic angular scatter of the 
cascade particles, which is estimated as Jiimit — "me/E,,. 
Comparing 5 from Eqs. ([5(1)1 and ([5^ with (Jumit one finds 
that the measurement of magnetic fields via the detection 
of cascade emission is possible for Bq larger than 



B, 



limit 



3 X 10-22 G 



{l + zf 



Ej 
0.1 GeV 



1/2 



A'b » 



1.8 X 10-20 G vm 

-1/2 







1/4 


[ A's 1 




[0.1 GeV J 




1 kpc 



X'r < De 



(41) 



^dclay, limit — 10 



gK(l 



(1 + ^) 



0.1 GeV 



-3/2 



(42) 



The time delay of the cascade induced tails of 7-ray 
flares depends only on the energy of the detected photons 
and on the distance to the source. This means in the case 
B <C Siimit one expects to 

• (a) detect echoes in the E^ ~ 0.1 — 1 GeV energy 
band in the flares of all blazars with spectra ex- 
tending up to the Eyg ^ 1 TeV energy band, 

• (b) detect the increase of the time delay of the 
echoes with the increase of the source redshift and 

• (c) detect the decrease of the time delay of the 
echoes with the increase of the 7-ray energy. 

Systematic observation of such delayed "echoes" in 
Fermi energy band would be an unambiguous evidence 
for the hypothesis of extremely weak or zero magnetic 
fields B ^ Siimit in the IGM. At the same time, non- 
detection of the 7-ray flare echoes with "universal" pa- 
rameters for all blazars would rule out the possibility of 
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FIG. 4: Left panel: grey shaded region shows the range of {B, As) parameter space which can be excluded via non-observation 
of "minimal possible" time delay of 7-ray flares by Fermi. Cyan colored limits on {B,Xb) are the limits shown in details in Figs. 
[T] and [21 Middle panel; region of (B, A_b) parameter space which can be probed via observations of time delays of 7-ray flares 
with Fermi (dark-grey shaded) and current and next generation atmospheric Cherenkov telescopes (orange hatched). Right 
panel: region of {B,Xb) parameter space which can be probed via observations of extended emission with Fermi (dark- grey 
shaded) and current and next generation atmospheric Cherenkov telescopes (orange hatched). 



extremely weak magnetic fields or, in other words, im- 
pose a lower bound on the EGMF strength. 



B> B 



limit ■ 



(43) 



and, in this way, rule out a range of " astrophysical origin" 
or " cosmological origin" models of the seed field. 

The curve B = Sumit for = 0.1 GeV (typical energy 
of 7-rays observed by Fermi) is shown by the thick black 
dashed line in the left panel of Fig. H] 



VII. MEASUREMENTS OF EGMF VIA TIME 
DELAYS IN 7-RAY FLARES 

A. Fermi 

li B > Biimit, the time delay of the 7-ray flare signal 
increases above the minimal possible value given by Eq. 

and starts to depend on B. Measuring the (energy 
dependent) time delay one can, therefore, determine the 
value of EGMF in the region of the size R ~ D^{E^^) 
around the 7-ray source. 

Typical durations of the fiares of blazars monitored by 
Fermi are Tflare ~ 1 — 10 days'^. If t{E^^,z) > 1, most 
of the primary source flux at the energy E.yg is absorbed 
in interactions with EBL and is subsequently re-emitted 
within the time Tdoiay at the energy E-y. The flux of the 



^ See the Fermi lightcurves of selected blazars at 
Ihttp : //f ermi . gsf c .nasa. gov/ssc/data/access/lat/msl_lc/| 



delayed emission could be estimated as 



-^delay 



Tdolay + Tfia 



Fo{{l + z)Ey 



(44) 



The spectra of all the TeV 7-ray loud blazars have pho- 
ton indexes harder than F = 2, which corresponds to the 
equal power emitted per decade of energy. This means 
that for the TeV blazars the energy flux Fq{{1 + z)Ej„) is 
normally higher than Fq{E^). This, in turn, implies that 
for sufficiently distant sources with t{E^^^ , z) > 1 (signif- 
icant fraction of the primary source power at the energy 
E^g is transferred to the secondary pairs) the delayed 
cascade flux is comparable or larger than the primary 
flare flux at the energy E^, as long as Tdeiay < Tflare- 
Thus, if the primary source flare is detectable with 
Fermi, the secondary cascade emission should be also 
readily detectable, at least for the range of time delays 

Tdeiay ~ Tflare ~ 1 — 10 d. 

To deconvolve the cascade contribution from the di- 
rect emission from the source one has to simultane- 
ously fit the lightcurves in different energy bands with 
a sum of direct time dependent emission from the 
source exp[—T{E^, z)] FQ{E^,t) plus a delayed contri- 
bution. The expected lightcurve of the delayed emis- 
sion is obtained via a convolution of the direct emission 
lightcurve with an exponential kernel [Tsj 



E 



delay(i,-B^) = / (l - fi-^^^^o Wq ( (1 + z)£;^„ , t') 

d!('45) 



exp 



t-t' 



Tciclay {E'y 



The GeV 7-rays are emitted by electrons/positrons 
produced by the primary 7-rays of the energies E-y„ ~ 
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1.3 TeV (see Eq. ((27)) ). Such 7-rays can travel over sev- 
eral hundred Mpc distances (see Eq. p^ ). This means 
that the timing observations with Fermi would reveal the 
magnetic fields averaged over many voids of the large 
scale structure. 

To estimate the maximal magnetic field detectable 
via the measurements of time delays, we notice that 
if the time delay becomes much larger than the dura- 
tion of the flare, the delayed flux diminishes by a factor 
Tdoiay/Tfiarc 3> 1, Compared to the primary source flux 
(see Eq. p5)) ). If the delayed flux becomes much smaller 
than the quiescent source flux, the delayed emission is 
difficult to detect. Assuming that the flux enhancement 
in the 1-10 d time scale flares of Fermi blazars is by a fac- 
tor of '-^ 10, we can estimate that the maximal detectable 
time delays are of the order of Tdoiay ^ 10^ days. Since 
the EGMF induced time delay is smallest at the largest 
energies, the maximal time delay is best detectable at the 
highest energies accessible for the observations, which are 
Ej ^ 100 GeV for Fermi. Substituting the estimated 
maximal Tdoiay and E-y into Eq. (j40|l one finds that the 
maximal magnetic field is 
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The range of EGMF strengths measurable via time delays 
with Fermi is shown as the lower dark grey shaded region 
in the middle panel of Fig. 01 



B. Cherenkov telescopes 

Ground based Cherenkov telescopes detect 7-rays of 
somewhat higher energies (Ej > 50 GeV) than Fermi. 
Next generation Cherenkov telescopes, like CTA, are 
expected to reach the low energy threshold of about 
E.y ~ 10 GeV. Higher threshold energy of Cherenkov 
telescopes leads to somewhat higher limiting magnetic 
field (SD 



It is clear that in the absence of magnetic fields in ex- 
cess of i3iiniit one expects to observe the "universal GeV 
tails of TeV flares" effect, discussed in the previous sec- 
tion also with the next generation Cherenkov telescopes. 
The prospects for detection and study of such " universal 
tails" with Cherenkov telescopes are potentially even bet- 
ter than with Fermi. The effective collection area of the 
ground based 7-ray telescopes at the low energy thresh- 
old is at least 4 orders of magnitude larger than the effec- 
tive area of Fermi. This enables detection of very short 
time delays. In fact, the measurements of time delay be- 
tween the signal of blazar 7-ray flares in different energy 
bands reach ~ 10 s already with the current generation 
instruments [6l,[6§|. At the energies above ~ 10 GeV the 
minimal measurable time delays with the next-generation 
instruments will be limited mostly by the intrinsic time 
spread of the 7-ray flare, rather then by the sensitivity 
of the instrument. At the moment, fastest observed vari- 
ability time scale of blazar flares is ~ 1 min, comparable 
to the minimal possible variability time scale, given by 
the light-crossing time of the blazar's central engine, the 
supcrmassive black hole [t^, [71|, . Assuming that the 
minimal measurable time delays are Tdoiay ~ 10 s, one 
can check that the fields of the strength at the level given 
by Eq. 0^1) could be readily measured. 

De Re-scaling the estimate of the maximal measurable 
magnetic field to the higher energies accessible for 
Cherenkov telescopes one finds that the timing mea- 
surements are sensitive to the magnetic fields with the 
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(47) 



and much shorter minimal delay time, Tiimit ~ 10^k(1 
rg-i)(l + z)-i [E^/IO GeV]"^/^ s. 



(48) 

The range of EGMF strengths measurable via time 
delays with Cherenkov telescopes is shown as the orange 
hatched region in the middle panel of Fig. [H 



VIII. MEASUREMENT OF EGMF VIA 
DETECTION OF EXTENDED EMISSION 
AROUND EXTRAGALACTIC POINT SOURCES 

A. Fermi 

If the magnetic field is large enough, deflections of the 
cascade e'^e~ pairs can be detected not only via time 
delay of the cascade emission, but also as extended emis- 
sion of the angular size ([55]) directly in the images of the 
extragalactic 7-ray sources. 

The extension of the source could be detected if the 
source size is larger than the point spread function (PSF) 
of a telescope. The minimal magnetic field measurable 
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via the study of extended emission is the one at which 
the extension of the source becomes larger than the tele- 
scope's PSF at low energy threshold. 

In the particular example of Fermi, the PSF depends 
on the photon energy, decreasing roughly as 6psF — 
0.3° GeV]"° ** below ~ 1 GeV and slowly im- 

proving from ~ 0.3° at 1 GeV to Opsf ^ 0.1° at 
~ 100 GeV. 

Substituting the reference energy E^ = 1 GeV and 
Gext = 0.3° into Eq 
measurable with Fermi 



The range of magnetic fields measurable via a study 
of extended emission around extragalactic point source 
with Fermi is shown as an upper dark shaded region in 
the right panel of Fig. 21 



B. Cherenkov telescopes 



^0 
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Angular resolution of Cherenkov telescopes is compa- 
one can find the minimal field rable to that of Fermi at the energies Ej ^ 100 GeV. 

Somewhat higher low energy threshold of Cherenkov tele- 
scopes results in a higher estimate of the minimal mag- 
netic field strength measurable via detection of extended 
emission, 
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(49) 

It is interesting to note that, in principle, the size of the 
extended source grows with the decrease of energy below 
E^ ^ 1 GeV, so that, apparently, it should be easier to 
detect the extension of the source at lower energies, close 
to the low-energy threshold of Fermi, E^ ~ 0.1 GeV. 
However, for Fermi, 0ps grows with the decrease of en- 
ergy below E^ ^ 1 GeV roughly in the same way as 
0ext, so that the decrease of the photon energy does not 
facilitates the detection of extended emission. 

The total fiux of the extended source does not de- 
pend on the magnetic field strength. This means that 
the surface brightness of the source decreases inversely 
proportional to the O^^^. Detection of extended emis- 
sion becomes impossible when the surface brightness of 
the extended source becomes comparable to the fluctu- 
ations of the diffuse background. If the minimal point 
source flux detectable by telescope is i^psF, the min- 
imal detectable extended source flux scales with the 
source size as Fext — (©ext/Qpsp)^- For Fermi, Fps 
depends on the photon energy approximately as Fps — 
10-12 (E^/l GeV) erg/cm^s at E^ ■> 1 GeV. 

The flux of the brightest TeV blazars reaches Fq 
10-1" erg/cm^s. If the source is at sufficiently large dis- 
tance, so that Tg{Ey^, z) > 1, the extended source fiux at 
the energy E^ is Fcxt ^ Fq . Such flux can be detected by 
Fermi up to the energies E^ ~ 100 GeV is the source size 
is not larger than Q^xt < 0.1° {E^/100 GeV)"\ This en- 
ables to estimate the maximal magnetic field measurable 
with Fermi as 



^0 



' 2 X 10-15(1 + 2)2 

Qcxt 



TU 



G, 



O.max 



o.r 

1.4 X 10-14(1 + z)i/2 



Abo 


-1/2 


©cxt 






1 kpc 




Lo.i° J 




.TO 



100 GeV 

G, 



3/4 



Ab > 



Ab < D, 
(50) 



2 X 10 

Qcxt 



i«r,(l + z)2 
G, 



10 GcV 



0.1 

2.5 X 10-1^9(1 



Ab > De 





-1/2 


©cxt 


1 kpc 




Lo.i" J 



1/2 



G, 



E~f 



10 GcV 



-1/4 



Ab < 
(51) 

Much larger effective collection area of Cherenkov tele- 
scopes explains their better performance at the highest 
energies. Assuming that photons of the energy up to 
^ 6 TeV (produced by the primary 7-rays of the energies 
E^g ~ 100 TeV could be detected, one can estimate the 
maximal EGMF strength from the condition that the size 
of the extended source should be smaller than the size of 
the telescope's field of view (FoV) 
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(52) 

Typical sizes of the FoV of Cherenkov telescopes are 
several degrees. Obviously, extension of the FoV up to 
~ 10°, expected e.g. with the high-energy part of CTA 
or with AGIS will extend the range of EGMF strengths 
accessible for the study toward stronger fields. 

The range of the field strengths measurable via study 
of extended emission around extragalactic sources with 
Cherenkov telescope is shown at the red right-hatched 
region in the right panel of Fig. HI 



IX. EVIDENCE FOR STRONG EGMF 

If the search of EGMF with {B, Ab) in the range out- 
lined by the grey shaded region in Fig. [2] will not give 
positive results and the lower bound on EGMF strength 
discussed in Section IVII would be confirmed, this would 
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signify that magnetic fields in excess of ~ 10~^^ G are 
present all over the IGM. A direct test of this result will 
be possible with the help of imaging of extragalactic 7- 
ray sources. Strong EGMF should cause a phenomenon 
of the presence of " universal" extended pair halos around 
all extragalactic multi-TeV 7-ray sources [73| , with prop- 
erties which depend on the source redshift, but do not 
depend on the EGMF strength. 

Namely, if the magnetic field is strong enough, tra- 
jectories of e"*'e~ pairs are strongly deflected over the 
cooling length De , so that the deflection angle S becomes 
S ^ 1. In this case the secondary IC photons are emit- 
ted isotropically from a region of the size R ~ D^{E^^) 
around the primary 7-ray source. If the primary source 
emits 7-rays isotropically, the expected flux of the ex- 
tended halo at the energy is equal to the primary 
unabsorbed source flux at the energy E^^. The angular 
size of the halo, 6'haio(£'.^) D^{E^^,)/D ~ 1/t0{E^„,z), 
depends on the source redshift, but not on the intrinsic 
properties of the source or on the EGMF strength. 

If the primary source (e.g. a blazar) emits the pri- 
mary 7-ray flux anisotropically, in a jet with an open- 
ing angle 6jotj the cascade emission fills a cone with an 
opening angle Oj^t and height H ~ (^70 ) i rather then 
forms a spherically symmetric halo around the source. 
In this case the size of extended source produced by 
the cascade emission is expected to be smaller, ^haio ~ 

Z?^(i?^Jejet/(i?-i?7(^7o)) -0jet/(Te(i?^o,^)-l)- At 

the same time, the fiux in the extended source is ex- 
pected to be suppressed by a factor 6?,^, compared to 
the primary source fiux at the energy E^^ . 

The halos around sources which are distant enough, 
T^E^g, z) ^ 1, and having significant fiuxes at high ener- 
gies, Fo{Ejg) ~ F{Ej), should be readily detectable with 
the current and next generation VHE 7-ray telescopes. 
In fact, if the primary source fiux is characterized by the 
photon index F < 2 the extended halo fiux could dom- 
inate over the primary source flux at the energy E-y, in 
the case of an isotropic primary source and be moderately 
below the point source flux in the case of an anisotropic 
primary source. For example, the flux from the brightest 
TeV blazars is at the level Fo{Ejg) ^ 10^^° erg/cm-^s at 
the energies E^g ~ 1 — 10 TeV. If the primary emission 
from the source is beamed into a jet with an opening 



angle Gjct ~ 5° (plausible assumption which implies the 
bulk Lorentz factor ~ 10 for the blazar jet), the expected 
EGMF-independent extended emission flux at the ener- 
gies i.e. at the level detectable by Fermi. 

The part of the {B,Xb) parameter space in which 
the existence of non-zero EGMF could be established 
via imaging 7-ray observations, but no measurement of 
EGMF is possible, is shown as an upper light-grey shaded 
region in Fig. [21 

X. CONCLUSIONS 

We have discussed the prospects of detection of weak 
magnetic fields in the intergalactic medium with the 
novel techniques of timing and imaging observations with 
ground and space-based 7-ray telescopes. These tech- 
niques enable measurements of extremely weak magnetic 
fields with strengths much lower than the ones accessible 
for the measurements with radio telescopes (via Zecman 
splitting and/or Faraday rotation techniques). 

We have demonstrated that using 7-ray observations 
one can detect, or rule out the possibility of existence of 
cosmologically or astrophysically produced "seed" mag- 
netic fields in the voids of the large scale structure, which 
arc conjectured to exist in a range of theories of the ori- 
gin of magnetic fields in galaxies and galaxy clusters (see 
Fig. 121). 

To summarize, we find that discovery or non-detection 
of weak magnetic fields in the voids of the large scale 
structure should provide, in the nearest future, a decisive 
test of the theories of the origin of cosmic magnetic fields. 
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